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Abstract. Starting with the Wigner distribution formulation for beam wave propagation in Holder 
continuous non-Gaussian random refractive index fields we show that the wave beam regime natu- 
(«~^ , rally leads to the white-noise scaling limit and converges to a Gaussian white-noise model which is 

' characterized by the martingale problem associated to a stochastic differential-integral equation of 

I the ltd type. In the simultaneous geometrical optics the convergence to the Gaussian white-noise 

model for the Liouville equation is also established if the ultraviolet cutoff or the Fresnel number 
' vanishes sufficiently slowly. The advantage of the Gaussian white-noise model is that its n-point 

correlation functions are governed by closed form equations. 



1. Introduction 



in 

m 
> 

' Laser beam propagation in the turbulent atmosphere is governed by the classical wave equation 

! with a randomly inhomogeneous refractive index field 

n{z,x) =n{l+n{z,x)), (z,x)g 

o ■ 



where n is the mean and n(x) is the fluctuation of the refractive index field. We seek the solution 
of the form E{t, z, x) = ^{z, x) exp [in{kz — cot)] + c.c. where E is the (scalar) electric field, k and 
w = kco/n are the carrier wavenumber and frequency, respectively, with cq being the wave speed in 
vacuum. Here and below z and x denote the variables in the longitudinal and transverse directions 
J3 ] of the wave beam, respectively. 

In the forward scattering approximation |^ , the modulation ^ is approximated by the solution 
^ I of the parabolic wave equation which after nondimensionalization with respect to some reference 

' lengths Lz and in the longitudinal and transverse directions, respectively, has this form 

■ (1) ik— + -A'^ + k^koL,h(zL„xL.^^ = 0, ^(0,x) = ^o(x) G L^fM'^), d = 2 

" " ' oz 2 

where k = k/ko is the normalized wavenumber with respect to the central wavenumber ko and 7 is 
the Presnel number 

= 

A widely used model for the fluctuating refractive index field n is a spatially homogeneous 
random field (usually assumed to be Gaussian) with the spatial structure function 

L'ndxl) =E[n(x + -)-n(-)]' = C'|x|2/3, \5l\ e {io, Lq), x = (z, x) G M'^+S d = 2 

where io and Lq are the inner and outer scales, respectively. Here and below E stands for ensemble 
average. 
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The refractive index structure function has a spectral representation 

sin (Ikl Ixl 



poo 

(2) Z?„(|x|) = 8^ / $„(|k| 

Jo 
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Ikllxl 



|kpd|k|, keM'^+^ 



with the Kolmogorov spectral density 

(3) ^nm\) = 0.033C2|k|-ii/3^ |k| e (^p^^p). 

Here the structure parameter depends in general on the temperature gradient on the scales 
larger than Lq. See, e.g., [22j, ^H] and [1] for more sophisticated models of turbulent refractive 
index fields. 

In this paper we will consider a general class of spectral density parametrized hj H €z (0, 1) and 
satisfying the upper bound 

(4) $(k) < i^(Lo2 + |k|2)--f^-l/2-«f/2 (l+^-2|k|2)-2^ k = k) G d = 2 

for some positive constant K < oo. Lq and £o in © are the infrared and ultraviolet cutoffs. The 
ultraviolet cutoff is physically due to dissipation on the small scales which normally results in a 
Gaussian decay factor [22j. We are particularly interested in the regime where the ratio Lq/Iq is 
large as in the high Reynolds number turbulent atmosphere. 

Let us introduce the non-dimensional parameters that are pertinent to our scaling: 

' Lx _ Lx _ Lx 

Lz^ -Lq ' 

In terms of the parameters and the power-law spectrum in @ we rewrite as 
az 2 ^ e 

with 

(6) /i 



(5) ifc^ + ^A^" + — ^F(-2,x)^'= = 0, ^'^(0,x) = ^'o(x) 



£3 



where a the standard variation of the homogeneous field n(2;, x) and V is the normalized refractive 
index field with a spectral density satisfying the upper bound 

(7) «>^,p(k) < K{r]^ + |k|2)-^^-i/2-rf/2 (1 + p-2|k|2) k G M'^+i, H e (0, 1) 

for some positive constant K. 

The generalized von Karman spectral density ^U], j22j 

(8) <i>,fc(k) = 2^-lr(F + l±i)^2/f^~(d+l)/2(^2 ^ |j^|2)-/f-l/2-d/2 

corresponds to the isotropic covariance function 

B,k{^) = E [y(x + ■)V{-)] = |??x|^E:j,(r?|x|), X = (z,x) G 
where -ftT/^ is a Bessel function of the third kind given by 



;>oo 

Kh{z) = / exp 
Jo 



dt. 



2 

For H = 1/2 we have the exponential covariance function B^k{5i) = exp[— r/|x|]. The additional 
ultraviolet cutoff imposed in the upper bound would then give rises to the covariance function 

S(x) = G*5,fc(x) 

where G is the inverse Fourier transform of the cutoffs. 
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For high Reynolds number one has Lq/Iq = p/r] ^ 1 and thus a wide range of scales in the 
power spectrum ((71). Note that in the worst case scenario the refractive index field loses spatial 
differentiability as p — > oo and homogeneity as rj ^ 0. The Gaussian field with its spectral density 
given by the right side of ((7| has H as the upper limit of the Holder exponent of the sample field. 
The Kolmogorov spectrum has the exponent H = 1/3. Since our result does not depend on d we 
hereafter take it to be any positive integer. 

Although we do not assume isotropic spectral densities, the spectral density always satisfies the 
basic symmetry: 

(9) ^(,,p)(e,k) = a>(,,,)(-e,k) = $(,,,)(^, -k), v(e,k) e R''+\ 

In other words, the spectral density is invariant under change of sign in any component of the 
argument because it is a characteristic function of a real-valued stationary process. 

We also assume that Vz(x) = V{z,x) is a square-integrable, z-stationary and x-homogeneous 
process with the (partial) spectral representation 

(10) y,(x) = j exp (ip • x)y,(dp) 

where the process Vz{dp) is the z— stationary orthogonal spectral measure satisfying 



(11) E y,(dp)y,(dq) =5(p + q: 



^{w, p)dw 



dpdq. 



We do not assume the Gaussian property but instead a quasi-Gaussian property (see Assumption 1, 
2 and 3 in Section 12.51 for precise statements) . 

If the observation scales Lz and are the longitudinal and transverse scales, respectively, of 
the wave beam then e ^ 1 corresponds to a long, narrow wave beam. The white-noise scaling 
then corresponds to e ^ with a fixed For convenience we set /x = 1. The white- noise scaling 
limit e — > of Eq. ^ is analyzed in (see also I^). The limit 7 — > corresponds to the 
geometrical optics limit. In this paper we study the higher moments behavior in both white-noise 
and geometrical optics limits by considering the Wigner transform of the modulation function. 

Our method is also suitable for the situation where deterministic large-scale inhomogeneities 
are present. One type of slowly varying, large-scale inhomogeneities is multiplicative and can be 
modeled by a bounded smooth deterministic function fj, = fj,{z, x) due to variability of any one of the 
three factors in © (see, e.g., |^, |21 for models with slowly varying a). The second type is additive 
and can be modeled by adding to fiV{z£~'^ a smooth background Vo(2,x). Altogether we 
can treat the random refractive index field of the general type 

, , ii(z,:x.) z , 
Voiz,x) + ^^^^V{-,x) 

e £^ 

with a bounded smooth deterministic modulation and background in the parabolic wave equation 
©. We describe the results in Section 12.31 but omit the details of the argument for simplicity 
of presentation. As the small-scale turbulent fluctuations are invariably embedded in a structure 
determined by large-scale geophysics this generalization is important for the practical application 
of the scaling limits. 

1.1. Wigner distribution and Wigner-Moyal equation. The Wigner transform of called 
the Wigner distribution, is defined as 

(12) Ty|(x,p) = ^ j e-P yvI/^(z,x + ^)^-(z,x-^)dy. 

One has the following bounds from (|12l) 

WIU < (27vr)-'^||^^(z,.)||L \m\\2 = {2^7:r"^W{zr)\\l 



|13j . |15| . [21j . The Wigner distribution has many important properties. For instance, it is real 
and its p-integral is the modulus square of the function cj), 

(13) / W^^(x,p)(ip = |M/^(x)|2, 

so we may think of Ty(x, p) as wave number-resolved mass density. Additionally, its x- integral is 

/ M^^(x,p)dx = (— )'^|$-|2(p/7). 

jRd 7 

The energy flux is expressed through (x, p) as 

(14) 1(^V^* -^'*V^') = / pT^=(x,p)dp 

2« jRd 

and its second moment in p is 

(15) y"|p|2VF(x,p)dp = |V^"(x)p. 

In view of these properties it is tempting to think of the Wigner distribution as a phase-space 
probability density, which is unfortunately not the case, since it is not everywhere non-negative. 
Nevertheless, the Wigner distribution is a useful tool for analyzing the evolution of wave energy in 
the phase space. Moreover, in the recent development of time reversal of waves in which a part of 
the waves is received, phase-conjugated and then back-propagated toward the source the refocused 
wave field is given by a Wigner distribution of mixed-state type (see (|25|) below) jjj, [22], [T2] . 

Moreover, the Wigner distribution, written as W^{x,p) = W'^{z,x,p), satisfies an evolution 
equation, called the Wigner-Moyal equation, 

(16) ^ + f ■ ^^^^^ + 7^'^^' = ° 

with the initial data 

(17) W^o(x, k) = ^ I e^'^ yM^ - ^)%ix + ^)dy , 
where the operator is formally given as 

(18) CIW! = iJe'^-^r'[W!{x,p + ^q/2)-W!{x,p-^q/2)]V{^,d<l) 

= 27-iy"wf(x,7q/2)Im [e-^27-^p-Xgiq-x-j>(^^^q) ^ 

We will use the following definition of the Fourier transform and inversion: 

•^-'^(x) = I e*P-"5(p)dp. 

When making a partial (inverse) Fourier transform on a phase-space function we will write J^i (resp. 
^i^) and (resp. to denote the (resp. inverse) transform w.r.t. x and p respectively. 

A useful way of analyzing C%W^ as formally given in (|18|) is to look at its partial inverse Fourier 
transform .F2~^£|VFf (x, y) acting on 

^2-iVFf(x,y)^ j e'P-yTy,^(x,p)dp = ^^(x + 7y/2)v&^*(x-7y/2) 



in the following completely local manner 

(19) ^2-i£^H^f(x,y) = -i7-i<^^F/(x,y)^2"iW^f(x,y) 
where 

(20) 57K'(x,y) = y/(x + 7y/2)-y/(x-7y/2) 

(21) n^(x) = y./,2(x). 

Hereby we define for every realization of the operator to act on a phase-space test function 
9 as 

(22) Lpi^, p) ^ -i7-'-^2 ['^7K'(x, y)-^2"'^(x, y)] 

with the difference operator 5^ given by H20|) for any test function 9 £ S where 

S = {0(x,p) G L\R^%:F^'eiK,y) G ^^(M'")} ■ 

We note that £| is skew-symmetric and real (i.e. mapping real-valued functions to real-valued 
functions). In this paper we consider the weak formulation of the Wigner-Moyal equation: To find 
G Z)([0,oo);L2(M2d)) ^^^^i that WW^h < ||VFo||2,Vz > 0, and 

(23) {W,, 9) - {Wo, 9) = ~k-' r {W^, p • Vx0) ds + - r {WIC19) ds. 

Jo ^ Jo 

Remark 1. Since Eq. \2'J\) is linear , the existence of weak solutions can he established straightfor- 
wardly by the weak-k compactness argument. Let us briefly comment on this. First, we introduce 
truncation N < oo 

Vn{z,^) = V{z,^), \V{z,^)\<N 
and zero otherwise. Clearly, for such bounded Vn the corresponding operator £| is a bounded 
self-adjoint operator on L^(R^'^). Hence the corresponding Wigner-Moyal equation preserves the 
Lp'-norm of the initial data and produces a sequence of -hounded weak solutions. Passing to the 
limit N ^ CO we obtain a Lp'-weak solution for the original Wigner-Moyal equation ifV is locally 
square-integrahle as is assumed here. However, due to the weak limiting procedure, there is no 
guarantee that the L'^-norm of the initial data is preserved in the limit. 

We will not address the uniqueness of solution for the Wigner-Moyal equation \2^) but we will 
show that as e — > any sequence of weak solutions to eq. i2cl\) converges in a suitable sense to the 
unique solution of a martingale problem (see Theorem 1 and 2). 

1.2. Liouville equation. In the geometric optics limit 7 ^ 0, if one takes the usual WKB-type 
initial condition 

*(0,x) = ^o(x)e^^W/^ 
then the Wigner distribution formally tends to the WKB-type distribution 

(24) Tyo(x,p) = |Ao|25(p-V5(x)) 

which satisfies J-2^Wq G L°°(M^'^). It has been shown 5 that the primitive WKB-type distribution 
H24|) can not arise from the geometrical optics limit (7 — > 0) from any pure state Wigner distribution 
as given by ()17|) but rather from a mixed state Wigner distribution of the form 

(25) T^o(x,k) = ^ J J e*k-y^o(x-^;«)^S(x+^;a)dydP(a), 

where P{ct) is a probability distribution of a family of states ^'q parametrized by a. The mixed 
state Wigner distributions generally give rise to a smeared initial condition, i.e. VFo(x, p) G L^(M^'^) 
even in the geometrical optics limit. This, instead of the WKB type, is the kind of initial conditions 
considered in this paper. 
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When acting on the test function space S, 11% as given by (|22|) has the following limit 

(26) lim r,0(x,p) = -T2 [Vxy.(x) • [^y^2"i0(x,y)]] = -Vxl4(x) • Vp0(x,p) 

in the L^-sense for all € 5 and all locally square-integrable Vz- Hence the Wigner-Moyal equation 
(|23|) formally becomes in the limit 7^0 the Liouville equation in the weak formulation 

(27) [Wl, 9) - {Wo, e) = k-^ r {W!, p • Vx0) ds-- r {W^,V^Vs • Vp0) ds, W G S. 

Jo ^ Jo 

The same weak-* compactness argument as described in Remark 1 establishes the existence 
of L^-weak solution of the Liouville equation except now that the operator (|26|) is unbounded 
and requires local square integr ability of VT4(-)- We will show that as e ^ any sequence of 
weak solutions of the Wigner-Moyal equation with any L^-initial condition converge as e, 7 ^ 
in a suitable sense to the unique solution of a martingale problem associated with the Gaussian 
white- noise model of the Liouville equation (see Theorem 2). 

In addition to the limit e — > we shall also let p ^ cxd and r/ — > simultaneously. We first 
study the case p — > 00, but ij fixed, as e — > 0. This means that the Presnel length is comparable to 
the outer scale. Then we study the narrow beam regime t/ — > where the Fresnel length is in the 
middle of the inertial-convective subrange. 



2. Formulation and main results 



2.1. Martingale formulation. The tightness result (see below) implies that for initial data 
the limiting measure P is supported in -L^([0, zq]; L^(M^'^)). For tightness as well as identification of 
the limit, the following infinitesimal operator will play an important role. Let = V{z/e'^, •) 
and zo < CO be any positive number. Let be the a-algebras generated by {Vf, s < t} and E| 
the corresponding conditional expectation w.r.t. J-'§. Let be the space of measurable function 
adapted to {J^z^z G M} such that sup^^^^El/zl < 00. We say fz G 'D{A^), the domain of A^ , and 
A'fz = gz if fz,9z^ and for /^^ = 6~^[Elfz+s - fz] we have 



supE|/f| < 

z,S>0 



limEl/f 

s^o ' 



9z 



00 



0, Vt. 



Consider a special class of admissible functions fz = f{{W^,9)),f'^ 
have the following expression from ()23() and the chain rule 



f'i{W!,0)),yf gc- 



we 



(28) A'fz = fz 

A main property of A^ is that 



Uw!,p-V.9) + ^{W!,C19) 
k £ 



(29) fz- A^'fsds is a J^l-martingale, \/f£V{A^). 
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Also, 
(30) 



Kfz - fs 



ElA'frdr ys<z a.s. 



(see C^). Note that the process is not Markovian and A^ is not its generator. We denote by 
A the infinitesimal operator corresponding to the unsealed process Vz(-) = V(z, •). 
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2.2. The white-noise models. Now we formulate the solutions for the Gaussian white-noise 
model as the solutions to the corresponding martingale problem: Find the law of Wz on the 
subspace of D{[0, co); L^(M^'^)) whose elements have the initial condition Wol^, p) S such 
that 

fi{Wz,e))- r \f'i{Ws,e)) \l{Ws,p-v^e) + P{Ws,Qoe 



is a martingale for each f £ C' 



+ Pf"{{Ws,9)){Ws,JCeWs) }ds 



with 
(31) 



JCeWs = J Q(0®0)(x,p,y,q)W^,(y,q)dydq. 



Here, in the case of the white-noise model for the Wigner-Moyal equation (Theorem 1), the covari- 
ance operators Q, Qq are defined as 



(32) 



Qo0 



^Jr(q)7"' [-2^(x, P) + 0(x, p - 7q) + ^(x, P + 7q)] dq. 



(33) Q(e®0)(x,p,y,q) = J e'^'<^-y^'^>^{ci')r^O{^,p - ^q' /2) - e{^,p + jq! /2)] 

X [9{y, q - 7q72) - 0(y, q + 7q72)] dq' 
and, in the case of the white-noise model for the Liouville equation (Theorem 2), 



(34) 



Qo^(x,p) 



Ap0(x,p) / cD^(q)|q|^dq 



Vqe(y,q), 



(35) Q(e®0)(x,p,y,q) = Vp^(x,p) 

?7 > 0, p < oo 
with the spectral density $^(q) given by 

$2°(q) = lim $^(q) = lim «>„,p(0,q), r? > 0. 

' p— >00 ' p— >CXD 

Note that the operators Q and Qq are well-defined for any test function 9 £ S in the former case 
for any H £ {0,l),r] > or ij = 0, H £ (0, 1/2), and in the latter case for H £ (0, l),Ori < p < oo 
OT H £ (0,l/2),0 = 7]<p<oootH£ (1/2, 1),0 < r/ < p = oo. 

To see that (|ITT |) - (|33|) is square-integrable and well-defined for any L^(M^'^) -valued process Wz, 
we apply J- 2^ to (|3H) and obtain 



(36) .F2-i/CeVr,(x,u) = .F2-i^(x,u) / e^^'-('^-y)$-(q')7 



gi7q'-u/2 _ g-i7q'-u/2 



X [0(y,q-7q72) -^(y,q + 7q72)] Wz{y,q)dydqdc(. 
(27r)2'^.F2-i0(x,u) j T2'9{y,y) [^-^VF,(y, y') - .^2"^VF.(y, -y') 



e-V-q'eV-(x-y)^oo(q/)^- 



gi7q'-u/2 _ g-i7q'-u/2 



dq'dydy'. 



The integral on the right side of is bounded over compact sets of (x, u) because 9 £ S, 
Wz £ L^fM^d) and the function 



is integrable in q' € M'^ and the associated integral is bounded over compact sets of u for any 
H £ (0,l),r7 > or rj = 0,H < 1/2. Hence the function on the right side of has a compact 
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support and is square-integrable. Similarly, one can show that H32() -H35 |) is well defined for H G 
(0, l),p <oo or H > 1/2, p = oo. 

In view of the martingale problem the white-noise model is an infinite-dimensional Markov 
process with the generator given by 



Afz = f's 



k 

This Markov process Wz can also be formulated as solutions to the Ito's equation 



(37) dWz = l^-^p • Vx + /c'QoJ Wz dz + kdBzWz, VFo(x) G L' 

or as the Stratonovich's equation 

dWz = ^T>-V^ + kdBzoWz, Wo(x) G l2(m2^) 
k 

where is the operator-valued Brownian motion with the covariance operator Q, i.e. 

E[dBze{x,p)dBz'e{y,q)] = 6{z-z')Q{e0e)ix,p,y,(i)dzdz'. 

Eq. dnZI should be solved in the space D{[0, oo); LI{R'^'^)), namely, to find G D{[0, oo); (1^2^)) 
such that for ah 9 G L2(IR2rf) 

(38) d{Wz,e) = (wz,(^lp-v^ + k^Qo^e^ dz + k{Wz,dBze), Wo{^) e l\r^''). 

Our results show that the solution to H38|) exists, is unique and satisfies the L^.i^ound 

(cf. Theorem 1, 2, Remark 1, 3 and Section f2.4|) . 

In view of H33|). (|32() . (|34() and (|35() we can interpret the white-noise limit e ^ as giving rise to 
a white-noise-in- 2; potential V* whose spectral density is bounded from above by 

i^*(r/2 + |k|2)-^*-'^/2 

for some constant K* < 00 with the effective Holder exponent H^, = H + 1/2 by observing that 

(39) lnn£|e(x,p) = -iJ^2[l-^5^V:{-K,y)J^^H{-K,y)\ , ^9 e S 

(40) hm £|e(x,p) = VxK*(x)-Vp0(x,p), V0 G 5 

£,7^0 

in the mean square sense. 

2.3. White-noise models with large-scale inhomogeneities. First we consider the case of 
deterministic, large-scale inhomogeneities of a multiplicative type which has fj,, given by as a 
bounded smooth function = ij,{z,x.). The resulting limiting process can be described analogously 
as above except with the term replaced by 

$-(k) ;u(z,x)^(z,y)$-(k), in Q 

$^(k) /z2(z,x)c^-(k), in Qo. 

As a consequence the operator Qo is no longer of convolution type. 

Next we add a slowly varying smooth deterministic background V(){z,x) to the rapidly fluctuating 
field fi{z , 'KjV {e~'^ z , 'k) . Namely we have 

, . /i(2;,x)^_ z , 

e 

as the potential term in the parabolic wave equation 



The resulting martingale problem has an additional term 

(41) - / k{Ws,CQe) ds 

Jo 

in the martingale formulation where CqO has the form 

£o^(x, P) = if e'^^'T"' [^(x, P + 7q/2) - ^(x, p - 7q/2)] Mz, dq) 

(42) = -il-':F2 [(Fo(x + 7y/2) - ^o(x - 7y/2))^2"^e(x, y)] 
for 7 > fixed in the limit, and the form 

(43) CoOix, p) = -VxFo(z, x) • Vp0(x, p) 
in the case of 7 — > 0. 

2.4. Multiple-point correlation functions of the limiting model. The martingale solutions 
of the limiting models are uniquely determined by their n-point correlation functions which satisfy 
a closed set of evolution equations. 

Using the function /(r) = r" in the martingale formulation and taking expectation, we arrive 
after some algebra the following equation 

(44) — — = rY^Pr ^-.^^"^ + E Qo(x,, p,)F(") + P 5] Q(x„ p„ Xfc, p,)F(«) 

for the n— point correlation function 

f(") (z, xi , Pi , . . . , x„, p„) = E [W^ (xi , pi) • • • W^{yin, Pn)] 

where Qo(xj,Pj) is the operator Qq acting on the variables (xj,pj) and Q(xj, pj, x^, p^) is the 
operator Q acting on the variables (xj, pj, x^, p^), namely 

Q(xj,Pj,Xfc,pfc)F(")(xi,pi) 



E 



I e^'i' (-y)$(,,^)(0,p)7-2 



n ^(x*,p* 

X [PF(x„ pj - 7q/2) - W{jcj,pj + 7q/2)][H^(xfc, pfc - 7q/2) - W{^k, Pk + 7q/2)] dq} • 
Eq. (|44() can be more conveniently written as 

(45) ^ = ^^p,-Vx^.f(") + P Q(x„p„Xfc,pfc)F(") 
with the identification Q(xj, pj, Xj, p^) = Qo(xj,Pj). The operator 

n 

(46) 2( 

Xj, Pj, Xfc, Pfcj 

j,fc=i 

is a non-positive symmetric operator. We note that the mean Wigner distribution can be exactly 
solved for from Eq. (|45|) for n = 1 and has a number of interesting applications in optics 
including time reversal. The 2-nd moment equation n = 2 is related to the problem of scintillation 
m (see, e.g., (H). 

The uniqueness for eq. 1)44(1 with any initial data 

fW(z = 0, XI, pi, . . . , x„, P„) = E [Tyo(xi, pi) • • • VFo(xn, Pn)] , G ^'(M'^) 

9 



in the case of the Wigner-Moyal equation can be easily estabhshed by observing that the operator 
given by H46() is self-adjoint. In the case of the Liouville equation, eq. (|45() can be more explicitly 
written as the advection-diffusion equation on the phase space 

f\Tp[n) 1 t2 " 

(47) = ^5:p,-V.,.fW + -5;I)(x,-x,):Vp,,Vp,fW 

j=l j,k=l 

with 



D(x,-Xfc) = y"e^'i'-(''^-"'=)$^(q')q'®qW 
DiO) = I <I>$;(q')|q?dq' 



with r/ > where D{0) is the Stratonovich correction term. In the worst case scenario the diffusion 
coefficient D{0) diverges as p — > oo but always well-defined as ?? ^ for H < 1/2. Moreover the 
diffusion operator 

n 

J2 D(xj-Xfc):Vp^,Vp, 

j,k=i 

is an essentially self-adjoint positive operator on C^(M^"'^) due to the sub-Lipschitz growth of the 
square-root of D(xfc — x^) at large |xj|, |xfc| jSI- 

2.5. Assumptions and properties of the refractive index field. As mentioned in the intro- 
duction, we assume that Vz{x) is a square-integrable, z-stationary, x-homogeneous process with a 
spectral density satisfying the upper bound 0. 

Let p{t) be a non-negative (random or deterministic) function such that 

|E[E,[y,(x)]E,[14(y)]]| = \E[E,[Vs{^)]Vt{y)]\ 

(48) < r{s- z)r{t- z)E[Vt^] , Vs,t > z,Vx,y G M"'. 

An obvious candidate for r(t) is the correlation coefficient defined as follows. Let J-'z and J-'^ 
be the sigma-algebras generated by {Vs : Vs < z} and {Vg : Vs > z}, respectively. The correlation 
coefficient r^^p{t) is given by 

(49) r^Ai)= sup sup E[hg]. 

E[hl=0,E[h2l = i " z+t 

E[9l=0,E[92] = l 

Lemma 1. The correlation coefficient rrj^p{t) as given by satisfies the inequality i4^ - 
Proof. Let 

/i,(x) =E,[y,(x)], gt{^) = Vt{^). 

Clearly 

gt G eL\P,n,J^+) 
and their second moments are uniformly bounded in x since 

E[/i2](x) < E[g!]{^) 

Eb'](x) = J m^)d^dci 

From the definition (|49|) we have 

\E[hs{^)htiy)]\ =\E[hsgt]\ < r,,pit - z)E^/^h',i^)]E^/^g^] . 
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Hence by setting s = t first and the Cauchy-Schwartz inequality we have 
E[/i2(x)] < rl^is - z)E[g^] 
E[hs{x)htiy)] < rn^p{t - z)rr^^p{s - z)E[g^], Vs, t > z, Vx, y. 

We assume 

Assumption 1. (i) The function r(t) in satisfies 

/ / E[r{s)r{t)]dsdt < oo. 
Jo Jo 

(ii) 

Corollary 1. Assumption 1 implies that the formula 

(50) y,(x) = / Eo [V;(x)] ds 

Jo 

defines a square-integrable x-homogeneous process. 
Proof. We only need to verify its second moment is finite: 



□ 



E 



= E 
= E 

< E 



oo roo 



z J z 
oo />oo 



E,[Vs{x)]E,[Vt{-K)]dsdt 
E,[Vs{x)]Vt{x)dsdt 



z J z 
oo fOO 



JO 



r{s)r{t)dsdt 



nvz 



which is finite by Assumption 1. 



□ 



One can adopt other alternative mixing coefficients to get the above results and higher order 
moment estimates, see Appendix A. Hereafter we will mainly focus on the correlation coefficient as 
it is most convenient to work with in the Gaussian case and we shall write explicitly the dependence 
of the correlation coefficient on r/,p as rrj^p{t). 

In the Gaussian case the correlation coefficient rri^p{t) equals the linear correlation coefficient 
given by 



(51) 
where 



rrj,p{^) = sup / R{t - Tl - T2,'k)gi{Tl,'k)g2iT2,'k)dkdTldT2 
91,92 J 

R{t,k) = j e**«<I>(^,,)(C,k)de 



and the supremum is taken over all (71,92 £ L?{W^~^^) which are supported on {—00, 0] X M"^ and 
satisfy the constraint 

(52) j R{t -t','k)gi{t,k)gi{t','k)dtdt'dk = j R{t - t' ,k)g2it,k)g2{t' ,k)dtdt'dk = 1. 
Alternatively, by the Paley- Wiener theorem we can write 

(53) r,Jt) = sup / e'«Vi(e,k)/2(e,k)$(e,k)dedk 

/l,/2 J 
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where /i, /2 are elements of the Hardy space Ti? of Li^iW^; <I>(,^ p))-valued analytic functions in the 
upper half ^-space satisfying the normalization condition 

j |/,(e,k)|2$(,,,)(e,k)dea!k = 1, j = 1,2. 

There are various criteria for the decay rate of the linear correlation coefficients, see |17| . 

Corollary 2. If Vz is a Gaussian random field and its linear correlation coefficient r^,p(t) is 
integrahle, then Vz is also Gaussian and hence possesses finite moments of all orders. 

This follows from the fact that the mapping from Vz to Vz is a. bounded linear operator on the 
Gaussian space. 

The main property of as a random function is that 

(54) AVz = -Vz, a.s. z£R. 

Since A commutes with the shift in x so the appearance of x in eq. is suppressed. 
We have the following simple relation 

1 



(55) 



lim E 

A— »oo 



Vzx{^)Vzx{y) 



lim 

A— >cx> 



3*(x-y)-p 



i^ 



1) ^{ri,p){^,P)d^dp 



Define the covariance functions 



and write 



i?,(x 



^ I e^(-y)-PcI.(^^p)(0,p)dp, Vz. 
Vz{^)Vz{y) 



E 



Bz{^) 



e^^-''^z{\^)dk 



where ^z{^) is its spectral density function. 

By the properties of the orthogonal projection Ez[-], we know that 



(56) 



E 



< E 



V{A)V{A) 



^(r,,p){i,^)didk 



Ez[V{A)]Ez[V{A)] 
for every Borel set A C M^^"^. 

Assumption 2. (i) The limit r^^oo = linip-»oo '^r^.p (*) exists and satisfies 



' ri,oo 



{t)dt 



lim 

P^'^ Jo 



r^^p{t)dt < oo. 



(a) As rj ^ 0, there exists a constant C such that 

rr,,p{t)dt < Ct]-^ I ri^p{t)dt 



For the Gaussian case with the von Karman spectrum (jHJ, a straightforward scaling argument 
shows that 

'rr,,p{t) = ri,p{rit) 
which is indepedent of p. This motivates Assumption 2. 
Set 

$-,(k) = |.,-.,(e,k) 
which is the spectral density of V^/(x) = F2/e2(x). 
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Define analogously to (|^ 
(57) ^|^(x, p) ^ -ir^T2 kF/(x, y)^2"'^(x, y) 



with 



<5^1//(x, y) ^ y/(x + 7y/2) - l7/(x - 7y/2). 



Lemma 2. (Appendix B) For each zq < oo there exists a positive constant C < oo such that 



sup E 

\z\<ZQ 

\y\<L 



(y) < C7^ min(7"\p) 



1 .i2-2H 



sup E 

\z\<zo 



sup E 

|2|<Z0 

|y|<i 



sup 

z|<zo 

|y|<-£' 



VyE 



y/(x) < Cri 

2 



-2-2H 



(y) < C'r/ ^7^1 min(/3,7 



-l\\2-2H 



(y) 



< Cr7-Vp^~'^|min(p,7~')| 



1m1-// 



sup E||p- Vx(/:",0)||^ < Cr^-'p 

\z\<zo 



-2 A-2H 



G S 



for all H G (0, l),e,7, < 1 < p, x, y € M"^ where the constant C depends only on zq, L and 9. 

We also need to know the first few moments the random fields involved. The case of Gaussian 
fields motivates the following assumption of the 6-th order quasi-Gaussian property. 

Assumption 3. 

(58) 

(59) 



sup E [5^y!{y)f < Ci sup E^ [S^V^f (y) 
\y\<L \y\<L 



sup E 

\y\<L 



(y) < C2 sup 

\y\<L 



(y) 



(60) sup E 

|y|<-f^ 



(y) < C3 sup E [5^V^f (y) sup 
\\y\<L / \|y|<L 



(y) 



for all L < CO where the constants Ci,C2 and C3 are independent of e,r], p,'j. 
From (|22() and (|57() we can form the iteration of operators >C|/2| 

£|>C|0(x,p) = -7-2^2 ky/(x,y)<5^y/(x,y)^2-ie(x,y 



The operator C.\C.%Q is well-defined if ^-^V^ and (5^V^/ are locally square- integr able. Other iterations 
of L% and £| allowed by Assumption 2 can be similarly constructed. 

The following estimates can be obtained from Lemma 2 and Assumption 3. 
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Corollary 3. 



E 



|£^,0(x,p)£^^(y,q) 



O I 7-* sup E|5^y/|2(y)E 
V lyl<^ 



E 



E 



I ne pe r\\\2 
\'-z*-'z^\\2 



O I 7-* sup E|5^l^/p(y)E 
V lyl<^ 



E 



O 7-^ sup E^ 

V |y|<i 

O ( 7"^ sup E^ 

V |y|<i 



(y)j =0(7?-Vm(p,7"')l'-''' 
(y)j =0(7?-Vm(p,7-')l'-''' 
0(^-Vm(p, 7-^)1'-''') 



■ E \b^V!f = O (r?-^| min (p, 7-')!'-'^) 



where the constants are independent of p, t], 7 and L is the radius of the ball containing the support 

ofT^^e. 

Assumption 4. For every 6 ^ S, there exists a random constant C5 such that 

(61) sup ||(5^V7J^2"^^ll4 < sup ¥}l'^\b-^V^{-K,y)\'^, G 5, e, ry, 7 < 1 < p 

z<zo ye zG[o,zo] 

|x|,|y|<i 

with C5 possessing finite moments and depending only on 9,zo, where L is the radius of the ball 
containing the support of J^2^6. cf. Lemma 2 and \6^) . 

For a Gaussian random field, Assumption 3 is readily satisfied for any a G (0, 1) by Lemma 2 
and a simple application of Borell's inequality ^ 

sup ||(5-^y/jr2""^0||4 < ||J^2" "^6*114 sup |(5-^V"/(x, y)| 

2G[0,zo] 
l=<:|,|y|<i 



z<zo 



(62) 



< C5log(3) sup EV2|5^v;-(x,y)'2 



e / z6[0,zo] 
|x|,|y|<i 



Vr/, 7 < 1 < /5 



where the random constants C5 has a Gaussian-like tail. 

Note that with 7 or p held fixed the first term on the right side of (|61|) is always 0(1). Compared 
to the corresponding condition ()62() for the Gaussian field condition H61|) allows for certain degree 
of intermittency in the refractive index field. 

Finally, we assume that for all p < cxd the refractive index field is smooth in the transverse 
coordinates. 

2.6. Main theorems. 

Theorem 1. Let be a z-stationary, x-homogeneous, almost surely smooth, locally bounded 
random process with the spectral density satisfying the bound Q) and Assumptions 1,2,3. Let 7 > 
be fixed. 

(i) Let r] be fixed and p be fixed or tend to 00 as e ^ such that 
(63) lim ep^~" = 0. 

Then the weak solution of the Wigner-Moyal equation with the initial condition Wq G 
L^(M^'^) converges in law in the space D([0, cxd); L^(M^'^)) of L"^ -valued right continuous 
processes with left limits endowed with the Skorohod topology to that of the corresponding 
Gaussian white-noise model with the covariance operators Q and Qo as given by Ii33\) and 
H&2\) . respectively (see also j| ) and \4^)- statement holds true for any H G (0, 1). 
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(ii) Suppose additionally that H < 1/2 and rj = 7]{e) such that 

(64) lim erj-^{r]-^ + p^'^) = 0. 

Then the same convergence holds. 

Here and below L^(R^^) is the space of square integrable functions on the phase space M^'^ 
endowed with the weak topology. 

The next theorem concerns a similar convergence to the solution of a Gaussian white-noise model 
for the Liouville equation. 

Theorem 2. Let he a z-stationary, x-homogeneous, almost surely smooth, locally bounded 
random process with the spectral density satisfying the bound Q) and Assumptions 1,2,3. 
Let 7 = 7(e) ^ as e ^ 0. Then under any of the following three sets of conditions 

(i) p < 00 and 77 > held fixed; 

(ii) H > 1/2, rj > fixed and p = p{e) — > cxd as e ^ such that 

(65) lim ep^-" = 0; 

(iii) H < 1/2, p < 00 fixed and rj = r]{e) such that 

(66) lim er/-2 = Q; 

£^0 

the weak solutions of the Wigner-Moyal equation ili6)] with the initial condition Wq G L'^{M?'^) 
converges in distribution in the space /^([0, 00); L^(M^'^)) to the martingale solution of the Liouville 
equation of the Gaussian white-noise model with the covariance operators Q and Qq as given by 
\cl4\ ) and ll&5\} . respectively (see also \41\ ) o.'i^d i4^)- 

Note that the Kolmogorov value = 1/3 is covered by the regimes of Theorem 1 and Theo- 
rem 2(i), (iii). 

Remark 2. As we have seen above, most of the assumptions here are motivated by the Gaussian 
case and we have formulated them in such a way as to allow a significant level of non-Gaussian 
fluctuation. In the Gaussian case, we have seen that a in i64\ ) OLnd \6b]) can be chosen to be any 
value in (0, 1/4). 

Remark 3. Both Theorem 1 and 2 can be viewed as a construction (and the convergence) of 
approximate solutions (via Remark 1) to the Gaussian white-noise models which are widely used in 
practical applications dS], j^. 

3. Proof of Theorem 1 and 2 
3.1. Tightness. In the sequel we will adopt the following notation 

(67) /,^/((Ty|,0)), f'^^f'{{wi,e)), f':^f"{{wi,e)), v/gC7-(m). 

Namely, the prime stands for the differentiation w.r.t. the original argument (not z) of /, /' etc. Let 
L denote the radius of the ball containing the support of T^^Q. Let all the constants c, c', ci, C2, . . . 
etc in the sequel be independent of p,rj,^ and e and depend only on zo,9., Il^^olb and /. 

First we note that since S is dense in L^(M^'^) and ||VF|^||2 < llW^olbjVz > 0, the tightness of 
the family of L^(M^'^)-valued processes {VF^,0 < e < 1} in D([0, 00); L^(M^'^) is equivalent to the 
tightness of the family in -D([0, co);S') as distribution-valued processes. According to jjlj, a family 
of processes {VF^,0 < e < 1} C D([0,oo);5') is tight if and only if for every test function 6 G S 
the family of processes {{W^, 9) ,0 < e < 1} C -D([0, 00); M) is tight. With this remark we can now 
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use the tightness criterion of [20j (Chap. 3, Theorem 4) for finite dimensional processes, namely, 
we will prove: Firstly, 

(68) lim limsupP{sup | (1^1,6') | > iV} = 0, Vzq < oo. 

N^oo z<zo 

Secondly, for each / G C°°(M) there is a sequence /| G 'D{A'^) such that for each zq < oo {A^ f^, < 
e<l,0<z<zo}is uniformly integrable and 

(69) limP{sup I/I - /((W|,^))| > (5} = 0, V(5 > 0. 



z<zo 



Then it follows that the laws of {{W^,6) ,0 < e < 1} are tight in the space of D([0,oo);]R) and 
hence {W^} is tight in L>([0, oo); L^(M2'i)). 

Condition (|H5|) is satisfied because the L^-norm is preserved. 

We shall construct a test function of the form /| = /z + /i ^ + /1 2 + /f 2- First we construct the 
first perturbation /f ^. Let 



Recall that 



Let 



(70) /!% = - / fAW!,Kcid)ds 



£ 



kef, (^T^'wi, ^-H^ MVs'^sJ'^^e 



2 zi 1 '-'•^1^ Z-' 2 

kefJwicie] 



be the 1-st perturbation of fz- 
Proposition 1. 



lim sup E|/i 2I =0, lim sup |/i 2! = in probability 

z<zo ' z<zo 



Proof. First 

(71) n\fiz\] < ewfiuiwohnmh 

< ce||/'||oo||H^o||2 sup [7-i5^y/(x,y) 
|x|,|y|<L L 

(72) = 0(er/-i|min(p,7-i)|i-^) 



which is of the following order of magnitude: 
(73) 



e, if ^) P held fixed 

e, if 7) ^ held fixed 

£r]~^ , if 7 or p held fixed 

e| min (p, 7^^)|^~^, if rj is held fixed. 
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and vanishes in the respective regimes. Secondly, we have 



sup \f!J 

z<zo 



z<zo 

< gup Ei/2|^-i<5^-l7/(x,y)|2 



(74) 



|x|,|y|<L 



by (|6ip . with a random constant c possessing finite moments, which is of the following order of 
magnitude: 

e^/^, if r], p held fixed 

if 7, 77 held fixed 
if 7 held fixed 
if p held fixed 
^ e^/2| min (p, 7^"'^)|"'^^^, if 77 is held fixed. 



(75) 



^1/2^-1 



which vanishes in the respective regimes. The right side of 1)741) now converges to zero in probability 
by a simple application of Chebyshev's inequality and assumption □ 

A straightforward calculation yields 



A'f! = -kefi(W!, 
where A'^O denotes 



k e 



A' 



Lie 



-/z{wi,ci 



+ kef:{w!,A'9)(w!,cie 



k 



-cie 



cf. (UHl)- Hence 

A' [fz + /: 



Iz 



= Al{z)+Al{z)+Al{z)+Al{z) 
where ^2 (-2) ^^d A^{z) are the coupling terms. 
Proposition 2. 

lim supE|A|(z)| = 



£^0 



z<zo 



Proof. By Lemma 2 we have 

l^ll < elirilooll^oll 
which is of the following order of magnitude 



\p-vM\^loh + \\p-v^{cie)\\2 



(76) 



ep 

ST]' 



2-H 
-1 



e{\mm{p,-f'^)\^-^ + p^~^] 



and vanishes in the respective regimes. 



if 1], p held fixed 
if ?7, 7 held fixed 
if p is held fixed 
if 77 held fixed 



□ 
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We introduce the next perturbations f2 fs z- 
(77) 4'^(0) = y"0(x,p)Qi(0®0)(x,p,y,q)</.(y,q)dx(ip(iydq 



(78) 

where 

(79) 

and 



Ql(0®^)(x,p,y,q) 



E 



Q'i0(x,p)=E ClCl9{^,p) 



where the operator £| is defined as in (|57() . Note that QiO and Q'i9 are 0(1) terms because of 
Clearly, we have 

A^^\cP)=K\{<p,Cie)(cP,Cl9] 



(80) 
Define 



2 fl/ 



E! 



ds. 



ds 



Let 
and 



Q2{e ^ ^)(x, p, y, q) = E [liei^, p)>C^e(y, q 
Q'2e(x,p)=E[£|r,0(x,p 



Let 




(81) 




(82) 




we then have 




(83) 


/l.z 


(84) 


f3,z 


Proposition 3. 





Q2^(x,p)(A(x,p) dxdp 



'^f\{wi,cicie)-4\w^,) 



Proof. We have the bounds 



lini sup E|/?^| = 0, lira sup \ f^^z\ =0) 3 — 2) 3. 



„ i "'?,zi 



supE|/2%i < supe^fc^iiriioo \\w4inc-ief2+m^?m)] 

z<zo z<zo L 



sup E|/|^| < sup e k 



2r.2ii i-/ 



"||P^o||2E||£^,£^,e||2 + E[4')(T^|)]" 

z<zo z<zo 

The first term can be estimated as in (|73() : the second term can by estimated as in ()73p by using 
the (jnU). 
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As for estimating sup^^^^ I//2I) J = 2, 3, we have 



sup|/|,J < supe^A^^liriloo \\Wof2\\Cie\\i + Af\wi 

z<zo z<zo 



27 2 n fi 



\\w,ucitM2 + Af\w, 



z<zo z<zo 

Using the assumption (|HT|) we can estimate the right side of the above as in ((7SJ. 



□ 



We have 



2,2 



P/;' -{W^,,ClO){W!,Cie) + A';'{W^,) +Rl{z) 



.(!)/ 



^'flz 



pfi -iwixim^A'^i^wi) +Ri{z 



with 



(85) 



k e 



-e'~k'f: 



(2) 

where Gn denotes the operator 



G 



(2). 



Q2{e ® 0)(x, p, y, q)(/.(y, q) dydq. 



Similarly 



(86) 

Proposition 4. 



k e 



k £ 



lim sup E|i?|(z)| = 0, lim sup E|i?|(z)| = 0. 

>^^0 Z<ZO e-^0 2<2Q 



WlClCie)-Af\w: 



Proof. Part of the argument is analogous to that given for Proposition |31 The additional estimates 
that we need to consider are the following. 
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In ii| (|H^: 



sup e E 

z<zo 



< ceV'l|W^o||2E{||VyVx^2"^e(x,y) 



X / E 



< ceV'llW^olbE 



5^y/(x,y)5^y/(x',y')J^2"^e(x',y')^2"^^|(x',y')dx'dy' 

2 



< ce^7 ^||ft'o||2 



VyVx^2^(x,y)E <5^y/(x,y) / |^2"^0(x',y')^2 T^|(x',y')| dx'dy' 

2 

2 



E||^2"'^-^2~'W^I||2 



< ce^l-^\\dUWo\\l 



+c\\eh\\Wo\\le'i-^ 



T^-^Vx- Vxe](x,y)E 
[^2"'Vx0](x,y)-VyE 



(y) 



(y) 



< cphWoWle'r^ sup E 
|y|<i 



< O min (p, 7-^)1'-'^ + min {p, T^')!'-^) 



(y)+c||0||2||H^o||ieV sup 

|y|<i' 

-i\a-H 



VyE 



(y) 



by Lemma 1^1 where L is the radius of the ball containing the support of 9. Further delineation 
yields 



sup e E 

z<zo 



g2^-2^1-/f 

e^p^"^! min (yO, 7~^)|"^~^ if ry held fixed. 



1], p held fixed 
if 7/, 7 held fixed 
if 7 held fixed 
if p held fixed 



Consider the next term: 



sup eE 

z<zo 



< ceV'l|W^o||2E{||5^y/(x,y)^2"'^(x,y) 



X / E 



< ceV'llW^olbE 



5^y/(x,y).F2-i0(x,y)E [5^V;^(x,y)] ' 

X I \j^^'e{^',y')^^'W^{^',y')\d^'dy' J 



< ceV'll^lbllTVoll^E 



6^V:{^,y)T^^eE 





2 







< 0(eV'|min(p,7-i)|3-3^^) 
by Corollary 13 



20 



In i?| 



sup eE 

z<zo 



< ell^olb sup -^/E 

z<zo 



(y)EV2|5^y/|2(y) 



= O ej-^ sup E 

= 0(er?-2|mm(/>,7~')l'-'''), 
by and Lemma [21 which can be further deUneated in the various regimes as 

e if p and tj held fixed 

e if 7 and r] held fixed 

er]~'^ if 7 or p held fixed 

e\m.m{p,"f^^)\^~^^ if r/ held fixed; 



(87) e^E\{W!,p-V^iQ'^e))\ < e^JE\{Wl,p ■V^iQ'S)\ 



< ceV'llW^oll 



Vy • VxE 



O { eS'2jE 



VyE 



<5^y/(x,y)J ^2-^e(x,y) 

2 



S.Vz 



(y) 



= 0(eVV-^|min(p, 7-^)1^-^) 

which in the various regimes has the following order of magnitude 

Hp and t] held fixed 
^2pi-H ^£ ^ ^j^^ ^ held fixed 

e^ri^"^ p^~^ if 7 held fixed 

e'^rj~'^ if /) held fixed 

e^/O"'^"''^! min (p, 7~"'^)|^^^ if rj held fixed; 



eE T^,^£^,Q'20 < eJ¥.\{Wl,ClQ, 



< ce^-=^||W^o||2E 



5^y/(x,y)E 5^y/(x,y) ^2"^0(x,y) 



O £^7"'^ sup E 

-2„-2| ^ -lM3-3H-j 



(y)Ei/2|5,y/|2(y) 



= O (e Tj I min (p, 7 
by Lemma 121 

Consider the test function fz=fz + fiz^f2z^f3z- We have 



1 w 



Set 

(90) R%z) = Rfiz) + R^2{z) + RU^), with i?^ (z) = A|(z) . 

It follows from Propositions El and jll that 

lim supE|i?^(z)| = 0. 



^.2 f aWi 



e^O 



z<zo 



□ 
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For the tightness it remains to show 
Proposition 5. {A^f^} are uniformly integrable. 

Proof. We shall prove that each term in the expression H89|) is uniformly integrable. We only need to 
be concerned with terms in R^{z) since other terms are obviously uniformly integrable because 
is uniformly bounded in the square norm. But since the previous estimates establish the uniform 
boundedness of the second moments of the corresponding terms, the uniform integrability of the 
terms follow. 

□ 

3.2. Identification of the limit. Our strategy is to show directly that in passing to the weak 
limit the limiting process solves the martingale problem formulated in Section 2.1. The uniqueness 
of the martingale solution mentioned in Section r2.4l then identifies the limiting process as the unique 
L^(]R^°')-valued solution to the initial value problem of the stochastic PDE (|37)1. 
Recall that for any C^-function / 



(91) 



Mm 



ft - / A' fids 
Jo 



fz + fl{z) + fliz) + fliz) - r jf, {W!,p- V, 

Jo k 



ds 







f':A^^\w!) + fiA\'\W^ 



ds 







R'{s)ds 



is a martingale. The martingale property implies that for any finite sequence < zi < Z2 < < 
■ ■■ < Zn < z, C^-function / and bounded continuous function h with compact support, we have 



(92) 



Vs > 0, Zi < Z2 < • • • < Zn < Z. 



0, 



•Afz = f's 



+ Pf'JA2{W, 



A2{e) = lim A'^\e) 



e), Ai{9) = lim A^^\9) = Qm 



Let 
where 

(93) . . . . , _ , , 

p — >00 p — fij^ 

as given in and respectively. For p— >oo,7— >Oase— >0 the limits in are not well- 
defined unless H £ (0,1/2) in the worst case scenario allowed by (O. Likewise, the convergence 
does not hold for G [1/2, 1) when r/ — > in the worst case scenario allowed by (jT)). 

For each possible limit process in Z)([0, oo); -L^(R^'^)) there is at most a countable set of discon- 
tinuous points with a positive probability and we consider all the finite set {zi, ...jZn} in (|^^ to 
be outside of the set of discontinuity. 

In view of the results of Propositions ^ |2 121 El we see that /| and ^^/| in (PTj) can be replaced 
by fz and Afz, respectively, modulo an error that vanishes as e — > 0. With this and the tightness 
of {W^f} we can pass to the limit e — > in (|92j) . We see that the limiting process satisfies the 
martingale property that 

, {Wz„,e)) [Mz+sio) - Mz{e)]} = 0, > o. 



E{h{{Wz„e),{Wz„ 

where 
(94) 

Then it follows that 



M,{e) = fz 



f 

Jo 



Afs ds. 



E[Mz+si9)-Mz 



Wu,U < z] 
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0, Vz,s>0 



which proves that Mz{6) is a martingale. 
Note that {W^,6) is uniformly bounded: 

\{wie)\<\\w42\\B\\2 

so we have the convergence of the second moment 

limE{(P7f,e)2}=E{(T^„0)2}. 

Using /(r) = r and in ()94|) we see that 

k 







ds 



is a martingale with the quadratic variation 



mW(0),mW(0) =P/ A2{Ws)ds = k^ {Ws,}CgWs)ds 
^ ^ Jo Jo 

where fCg is defined as in (|ST|) . 

Appendix A. Mixing coefficients and moment estimates for Vz 

Let .7^2 and J^^ be the sigma-algebras generated by {Vs : Vs < z} and {T^ : Vs > z}, respectively. 
Consider the strong mixing coefficient 



a{t) = sup sup \P{AB) - P (A) P{B)\ 
1 



2 AeJ^, 



sup ¥.[\P{A\Tz) - P{A)\] 



z + t 



which can be used to bound the first order moment: 

(in], Corollary 2.4). Hence the integrability of a{t) implies that Vz has a finite first order moment. 

To bound the higher order moments of Vz one can consider, for example, the general L^-mixing 
coefficients 

</.p(t) = sup ¥}I^[\P{A\J^,)-P{A)\^], pG[l,oo) 

= sup sup E[%], p^'^ + q^^ = l, pG[f,oo) 

We note that a{t) = (j)i{t) and for p = oo 

(pooit) = sup sup |P(^|-B) - P(^)|, Vt>0 

^^•^t+z P(B)>0 

= sup ess-sup^|P(A|.F,)-P(^)| 



A^K+z 



is called the uniform mixing coefficient jH]. In terms of (j)p one has the following estimate 

(95) |E[/ii/i2] -E[/ii]E[/i2]| < 2°'''^('?,2)^^(^)i/«]Ei/M[/j^P]]Ei/g[/ja] 

for u,v,p,q £ [l,oo],u^^ + v^^ = l,p^^ + Q^^ = 1 and real-valued hi G L^{i^,J^z, P),h2 G 
L^P{Q,T^_^_^, P) (see (HI, Proposition 2.2). In particular, foi q > 2,v = q/p, 

(96) |E [hih2] - E[/ii]E[/i2]| < 4(/)p(t)('?-P)/mi/^[/i^]EV'?[/i?], ^j-^ + = 1 
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by which, along with the Holder inequality, we can bound the second moment of Vz as follows: 
First we observe that for s,t > z and hi = ¥,z{Vs), h2 = Vr 

E [E,[14(x)]E,[K(x)]] = E [E,[y,(x)]V;(x)] < A^r - z)('?-P)/mi/«[V7]Ei/'?[E^ [Vs]]. 
By setting s = t first and the Cauchy- Schwartz inequality we have 
E[e2[14]] < 4,/.p(s - z)(5-P)/'?E2/g[y^g] 
E[E,[l^,(x)]E,[V;(x)]] < 4,/.p(s-z)(^-*')/(2'?)^p(r-z)(''-P)/(2'?)E2/5[y/], s,t>z. 

Hence 

/"OO roo roc roc 

E[V^] < 2 / E[E,[Vr]K,[Vs]]dsdT + 2 / E [Ko[Vr]Eo [Vs]] dsdr 

Jz Jz Jo JO 

/ i-oo \ 
< 8E2/'?[y/] / </.p(t)(9-P)/(2'?)(it 



OO 

2 



< 8E'/^[Vz']y^°^<l>l%{t)dt 



2/5 

which is finite if 4>Q^^{t) is integrable (if Vz is assumed to have a finite 6-th order moment). 
When Vz is almost surely bounded, the preceding calculation with p = 1, q = oo becomes 

E[V^] < 8 hrn Ei/«[V7] 

1/2 

which is finite when (p-^ (t) is integrable. 

In order to bound higher order moments in the non-Gaussian case, one can assume the integra- 
bility of the uniform mixing coefficient (j){t) = (pooit)- Then we have 

\P{A\J^z)-PiA)\<cPis-z), yAeJ^s, s>z 

and for p G [1, oo),^^"'^ -|- q^^ = 1 

(97) |E [VslJ'z]] < 2^/P^^/P{s - z) [E ^!\Tz] + E [y/]]^/" 

(cf. [HI, Proposition 2.6). Using (|^7)) and the Holoder inequality repeatedly we obtain 



(98) E<^ / E'^s\Tz]ds\ < c 



IP 



(j){s)ds 



E[vn. 



Hence the integrability of (j){t) implies that Vz given by (|5U() has a finite moment of any order p < oo 
if Vz has a finite moment of order p. 
In summary we have 

Proposition 1. (i) Assume that E[Vz] < oo,p £ [l,oo). // the uniform (L°^-) mixing coeffi- 
cient (poo{t) ofVz is integrable then Vz has finite moments of order p. 

(ii) Assume that E[Vz] < oo- If the 2/ 5-power of the L^^^ -mixing coefficient (j)Q/^{t) is integrable, 
then Vz has finite second moment. 

(iii) Assume Vz is almost surely bounded. If the square-root of the alpha- (L^-) mixing coefficient 
(j)i{t) is integrable then Vz has finite second moment. 
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(i) Estimation of sup|z|<zo IE 
|y|<i' 



Appendix B. Proof of Lemma 2 
{SjV^'f] (y) : We have that for 7p < 1 



sup E 

|z|<20 



(5^y/(x,y))^ 



< 



< 



sup f 4| sin (7y • k/2)|2$(^^^)(^, k)d^dk 
\z\<zo J 

sup / |7y •k|2$(^^^)(^,k)dCdk 

|z|<zo J 

sup / (r;2 + |k|2 + |C|2)-^-('^+i)/2|k|'^+id^d|k| 
\z\<zoJ\k.\<p 

< CiotVI'sUP / / (772 + |k|2 + |^|2)-^^-('^+l)/2|k|'^+ld|k|dC 

< cinVl' sup / + |cn-^+V2rf^ 

UI<2n^l£l</9 



< cW\y\' [ \r'"^'d^ 

For /97 > 1 we divide the domain of integration into Iq = {|k| < 7~^} and Ii = {|k| > 7~^} and 
estimate their contributions separately. For Iq the upper bound is similar to the above, namely, we 
have 

sup / 4|sin(7yk/2)|2$(^,^)(e,k)dCdk < C137VI' 7"'+''') • 
For Ii we have instead that 



sup j 4|sin(7yk/2)|2$(^^^)(^,k)dedk 

|z|<20 J h 

< 4 sup / ^^r,,p){^,k)d^dk 

\z\<zo J h 

< CM sup/ / {v' + m' + \^\Y"-^'^'^/'\M'-'d\m 

\z\<zoJ\^\<pJ\\^\<p 

< sup / + ler^'-'^'d^ 

\z\<znJ\£\<p 



< C16 / 

Put together, the upper bound becomes 



2H-1 



sup E 

Z|<ZQ 

|y|<i 



(5^y/(x,y))2 < C7'|min(7-i,p)|' ^,v<l<p. 
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(ii) Estimation of sup|_2|<2,Q IE V^/(x) : It follows from the argument for Corollary 1 that 



E 



2H 



(y): First note that the correlation coefficient for SjV^ 



(iii) Estimation of sup|z|<zQ E 
|y|<i 

is bounded from above by cr{t) for some constant c > 0. Then we have as above that 



E 



ri,p{t)dt\ 7^ |min(7"\p)| 



(iv) Estimation of sup|z|<zo 
|y|<i 

ceding calculation we have 



VyE 



(y) 



: By the Cauchy- Schwartz inequality and the pre- 
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sup 


VyE 




(y) 


Z|<ZQ 







< C3 

< C4??" 



VxF=(x + 7y/2) + VxF^(x - 7y/2) 

\ 2 



ri,p{t)dt 



-i^p" ^|min(p,7 ^)| 











x,y)]' 




-H 





(v) Estimation of supj^j^^o 

E||p • Vx(>C|6')|||: A similar line of reasoning and a straightforward 

spectral calculation yield that 



E||p- Vx(/:|^)||2 = E||Vy • Vx7"M^y/J7^^||^ 
< ciE||VxV;"j^2"^e||i 



< 



< C3?7 



4-2// 



26 



Appendix C. Proof of Corollary El 
By the Cauchy-Schwartz inequality we have the following calculation: 



E 



E 



E 



+ E 



E 



E 



5^y/(x,x')<^^F/(y,y')J^2"'^(x,x')^2~'^(y,y': 

2^ 



e|2 



-1/ 



E 



-1^ 



and 



and 



O I sup E|(5^F/|^(y)E 
\\y\<L 



(y) 



E 



C[-i-^ \ I E E U^Vi] {J^^^e^dxdy + 



•^2''^(x,y] 



O sup E|(5^y/r(y)E 
Vlyl<-f' 



E 



< C2|7 
= C27-' 



-4 



E 



E 



E 



<^7K' 



E 



cicie{x,p) 



e) (ix(iy + 



E 



^2^e{x,y) 



\ sup E^ 




\\y\<L 





(y) 



where Ci,C[, C2 are constants independent of p, ry, 7 and L is the radius of the ball containing the 



support of ^2 Similarly we have that 



E 



O sup E^ 



\\y\<L 
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